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INTRODUCTION 

It is well known that the modulation function of an in- 
commensurate modulated phase is smooth at the transi- 
tion from disordered to incommensurate, and that it may 
become discontinuous further from the transition line. In 
the incommensurate phase it marks the transition from 
a simple modulation to a structure that is locally like the 
commensurate structure that usually appears at another 
phase transition from incommensurate to commensurate 
modulation. The phenomenon of the appearance of a dis- 
continuity is often called the transition by breaking of an- 
alyticity [1] or the discommensuration transition [2], the 
latter because at the transition the structure is locally 
periodic but with domain walls, called discommensura- 
tion. The transition is accompanied by an opening of a 
gap in the phason frequencies [3] . It has been studied on 
several models [2-9] , like the discrete frustrated (p'^-mode\ 
(DIFFOUR), the Frenkcl-Kontorova model and a double 
chain model for incommensurate composites. 

From these studies the qualitative change in the modu- 



lation and the phason gap opening have been well estab- 
lished. It is not so clear whether the transition is a tran- 
sition in the thermodynamic sense, and whether there 
are discontinuities in the thermodynamic quantities. To 
that end we analyze the situation for a displacively mod- 
ulated structure using the DIFFOUR model. First the 
ground state properties are analyzed. Then temperature 
effects are studied. The analytical results are checked 
with numerical calculations. 



DEFINITIONS AND GROUND-STATE 
PROPERTIES OF THE MODEL 



The model under consideration is the 3D DIFFOUR 
array. It is a cubic lattice of 2-4 anharmonic double 
well oscillators, with a harmonic coupling. Each atom 
is coupled with 6 nearest neighbours and 2 next-nearest 
neighbours in the z-direction. The potential energy of 
the model is expressed as 
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Here A, B, C, D are constants. The third sum (marked with a prime) runs over all pairs of the nearest neighbours 
n and m, and the last sum (with a double-prime) is over the pairs of the next- nearest neighbours in the z-direction. 

The ground-state properties of the model depend on the two dimensionless combinations: a = —A/C and d = D/C. 
Along a line in the a — d plane, the modulated phase with an incommensurate modulation period becomes energetically 
favourable: 



(2) 



where / is a periodic function with, generally, an irrational period. For a small a, /(z) is a sinus-like continuous 
function. With an increase of a the modulation function becomes less smooth and at certain temperature can loss 
the continuity ("breaking of analyticity" or discommensuration transition). The part of the ground-state "phase 
diagram", containing the incommensurate phase, is schematically shown in Fig.l. Typical profiles of the modulation 
function are shown in Fig. 2. The plateaus at the discontinuous curve correspond to a locally periodic structure with 
discommensuration. 
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One can note that even the discontinuous modulation function roughly resembles the graph of sinus. Therefore it is 
reasonable first to consider a sinusoidal trial modulation fo{z) = xq sin kz. The minimisation of the potential energy 
gives cosfc = xl = ^ - AC-8D). 

The following procedure allows to describe qualitatively the appearance of the discontinuities. Let us use fo{z) as 
a starting point, and find a function fi{z) that minimises the potential energy of each layer in the fixed field of its 
neighbours: 

Vi{z) = + + 4C + 4i?)^ - 2hiz){Cfo{z + l) + CMz-l) + Dfoiz + 2) + DJ^{z - 2)) = min . (3) 




FIG. 1: A part of the ground-state phase diagram for the 
DIFFOUR model (1), including the incommensurate (IC) 
phase. Regions denoted by N=l and N=6 correspond re- 
spectively to the (commensurate) ferroelectric phase and to 
the modulated phase with a integer period of 6 layers. 
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FIG. 2: Examples of the ground-state modulation function 
above (a = 0.3) and below (a = 1.0) the breaking of continuity 
transition. Both curves are plotted for d = —0.4. 



One can consider this expression as the first iteration 
in the determination of the modulation function for a 
fixed period of modulation. However, we shall consider 
only the first iteration. Physically, this is a condition 
for the ground-state state of an anharmonic oscillator 
in the given external field formed by the neighbouring 
particles. The value of the field passes through zero at 
z = and at z = n/k. The qualitative behaviour of 
/i (z) is determined by the sign of — ^ -I- 4C -I- 4D. If this 
combination is positive, fi{z) is a continuous function. 
Otherwise, discontinuous steps from —\/A — AC — AD to 
y/ A ~ AC — AD appear at z = and at z = 1/2. 

The situation at the points of discontinuity resembles 
quite much the Landau scenario of second-order phase 
transitions. Thus the contribution to the ground-state 
energy from these two points is of the Landau type - for 
example, it shows a discontinuity in the second derivative 
with respect to A. However, to obtain the total energy, 
one should integrate over the modulation period. It is 
useful to consider the derivative of the ground-state en- 



ergy: 



= [iMdz (A) 

A dA ^ J 2 ^ ' 



To establish this formula we took into account that 
= 0, as Vgs should be minimised with respect to 
Xm. The behaviour in the transition point with / = fi{z) 
can be analysed by means of the expansion in powers of 
{-A ■^AC + AD). It turns out that both dVgs/dA and 
d'^Vgs/dA^ are continuous, and only the third derivative 
shows a step. 

Quantitatively, the approach described above is very 
inaccurate. For example, for D = —Q.AC it predicts the 
breaking of continuity a,t A — 2.AC, whereas actually 
it takes place a.t A ^ 0.6C. One can increase the ac- 
curacy, considering instead of lj2Jl the minimum of each 
three neighbouring layers in a given field of the surround- 
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ing layers: 

V3^{ftiz) + fH^ + l) + ftiz-l))/4 + 

+ 4C + 4Z?)(/2(z) + + 1) + /2(z _ l))/2 
-2C/i(z)/i(z + 1) - 2C/i(z)/i(z - 1) 
-2D/i(z-l)/i(z + l) 

-2D/i(z)(/o(z + 2) + /o(z-2)) (5) 
-2A(z + l)(C/o(z + 2) + Dfoiz + 3)) 
-2/i(z - l)(Cfo{z - 2) + Z?/o(z - 3)) = min. 

The function fi becomes discontinuous when the minimal 
eigenvalue of the 3x3 matrix of the second derivatives 
of 1/3(0) becomes zero. This gives the following condi- 
tion for the breaking of continuity: —A + 4C + 3-D — 
+ {D - 3f'fil))'^) + 3fl{l)/2 = 0. The prediction for 
the D = — 0.4C system is now much better: A « l.OC. 

It should be noted also that the continuity of the mod- 
ulation function actually breaks down not at the single 
point z — 0. Indeed, the exactly calculated modulation 
function is a fixed point of the transformation Q. It is 
clear from this transformation, that if there is a discon- 
tinuity at certain z, the discontinuity at z ± 1 should be 
present also. Therefore, in the mathematical sense, f{z) 
is discontinuous in all points when the continuity is bro- 
ken in a single point. Practically (see Fig. 2), the "parent" 
discontinuity at z = 0,7r/fc is the most remarkable, the 
magnitude of "satellite" discontinuities falls down fast. 

Although the above approach is very simplified, it 
clearly shows the physical picture of the breaking of con- 
tinuity. The qualitative prediction about the continuity 
of (fVgs/dA^ seems to be correct as well. Fig. 3 shows a 
numerical result for this quantity, which is a smooth line. 

At the same time, it is worth to note that the "break- 
ing of continuity" affects the spectrum of elementary ex- 
citations essentially. For a continuous modulation func- 
tion, a phason branch is present in the spectrum [2, 3, 9]. 
It shows a linear (phonon-like) dispersion, with a phase 

velocity w (X (^^^{df/dz)^dz^ . The zero-frequency 

Goldstone mode corresponds to the uniform shift of the 
modulation phase; its eigenfunction is uo cx vdf /dz (the 
normalisation is J v?dz = 1). For a discontinuous f{z), 
the integral for v diverges, indicating the presence of a 
gap at the beginning of the spectrum. This regime can 
be described as follows. Formally, one can still calcu- 



late uq from the above formula. For a modulation func- 
tion one obtains that uq becomes the two singular 
peaks at 2 = 0, 7r/fc. Physically, this means that the 
low-lying excitations are the vibrations of atoms only in 
the points where the modulation is broken. Therefore, 
their frequency is determined by the second derivative of 
Vi, cx d'^Vi/dz'^. One can interpret these vibrations 
as domain wall oscillations. Consequently, the presented 
model predicts a gap, which is proportional to the square 
of the discontinuity of the modulation. 
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FIG. 3: The second derivative of the ground-state potential 
energy with respect to a for the system with d — —0.4. There 
is a step at a = —0.45, corresponding to the transition from 
para-phase to the incommensurate state. No singularities at 
the breaking of continuity point (a « 0.6) are observed. 



FINITE TEMPERATURE: THE MEAN-FIELD 
APPROACH AND MONTE CARLO RESULTS 

Let us turn to the case of finite temperature T. The 
partition function is equal to 

Z = J exp{-V/T)dxi...dxn (6) 

(Boltzmann constant is equal to unity). 

In the mean-field consideration, the inter-atomic cor- 
relations are neglected, so that the average position of an 
on-site oscillator x{z) is given by a self-consistent equa- 
tion 



x{z) = x{f{z)) 



dF{f{z)) 
Of 



(7) 



-Tin / cxp ( -i 



{-A+12C + 4D)x'^ 



F{f{z)) 

}{z) = 2C{'ix{z) + x{z-l)+x{z + l)) + 2D{x{z~2)+x{z + 2)). 



f{z)x dx; 
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Here f{z) is an average force acting on the oscillator from its neighbours. 

Note that these equations determine the one-to-one correspondence between x{z) and f{z), so that one can construct 
a function f{x). It can be checked that equations iQ are the minimum conditions for the function 

Fmf = J2 i^ifM - F{0) + Xnf{xn) - (6C + 2D)xl)) + (8) 

n 

/ // 
C ^ ^ {Xn •^ra) ~t- D ^ ^ {Xji Xjn ) . 

Thus in the mean-field approximation the problem at a finite temperature is replaced by the problem for the ground- 
state on an effective system (0. The on-site potential of this effective system depends on temperature. The value 
of Fm f can be interpreted as a mean-field expression for the free energy of the system. As the temperature goes to 
zero, ® tends to (jTJ, if A is not too large {A < 12C -I- AD). 

The breaking of continuity of the modulation function occurs at a certain critical temperature T^. The considerations 
from the previous section can directly be applied here. In particular, the heat capacity, being the second derivative 
of the free energy, does not show a discontinuity at Tc. 

To go beyond the mean-field approach, it is necessary to consider the role of the fluctuations. Here we discuss, if 
the fluctuation phenomena can destroy the phase with a broken continuity. As it is demonstrated, the breaking of 
continuity occurs due to the appearance of the effective double- well potential for the layer with z — (formula Q)- 
To incorporate the fluctuation in this model, it is natural to suppose that in the zeroth approximation the modulation 
function is sinusoidal and fluctuations are absent. Instead of © consider the statistical sub-ensemble with energy 

E X + + E ^(^' - ^j)' - 2^ ^^(^/o(^ + 1) + Ck{z - 1) + DUz + 2) + DUz - 2)) (9) 

i i,j i 



Here the indices run over the layer with a given z- 
coordinate. The layer is in fact the 2D discrete (/)^ model 
[10] placed in an external field. The dependence of the 
average of Xi on the value of the external field is discon- 
tinuous, if the layer without the external field shows a 
phase transition. Since a phase transition occurs in the 
2D discrete (j)^ model indeed, we conclude that fluctua- 
tions do not destroy the discontinuity of the modulation 
function for the model under consideration. It can be 
noted that for a 2D DIFFOUR model this is not true. In 
fact, in this case we would come to the ID discrete 0^ 
chain, which does not show a phase transition. 

To check the conclusion about the thermal behaviour 
of the heat capacity, we performed a Monte Carlo simu- 
lation of the system with a modulated ground state with 
a discontinuous modulation function. The system was 
studied for a = 1.5,6? = -0.4. A slab of 8 x 8 x 557 
atoms with periodic boundary conditions and the stan- 
dard Metropolis algorithm was used. Initially we put the 
system in its ground-state, and set the temperature to be 
small. Then the temperature was slowly increased (about 
3 • 10^ Monte Carlo trial steps were used for each value of 
temperature). The average potential energy < E > was 
calculated at each point. Afterwards, the thermal capac- 
ity d < E > I dT was obtained. Thermal behaviour of 
this quantity is presented in Fig. 4. There is a clear peak 
at TB /C^ « 3.7 corresponding to the transition between 
the modulated and disordered phase. The breaking of 




TB/C" 

FIG. 4: Monte Carlo result for the heat capacity of the system 
a — 1.5, d — —0.4. The ground state of this system is incom- 
mensurate with a broken continuity. The data shows a clear 
peak at the transition to the disordered phase {TB/C^ ~ 3.7) 
and no other pronounced features. 



continuity of the modulation function occurs at a lower 
temperature, but it does not result in any remarkable 
feature oi d{E{T)) / dT . 
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CONCLUDING REMARKS 

In several model systems for aperiodic crystals there 
is a transition from a smooth modulation function to a 
modulation function with discontinuities. This transition 
has been called the discommensuration transition of the 
transition by breaking of analyticity or continuity. In the 
phase space region where the modulation function is con- 
tinuous there is in these models a phason excitation of 
zero frequency. In the other region there is a phason gap. 
However, neither the structural transition nor the open- 
ing of a phason gap has been observed in experiment. 
There are even few systems where a gap-less phason or 
phason with very low frequency has been found. Gen- 
erally this is attributed to pinning by defects, or other 
friction mechanisms [11]. Another possible method to 
check the existence of the transition is by measuring ther- 
modynamic quantities. In the model studied here, the 
DIFFOUR model, the transition is very weak. It can be 
called to be of third order. The conclusion is that the 
transition could only be observed by a precise structure 
determination. However, in general the discontinuity will 
be not strong enough to be seen. 
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